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SPECTRUM AND EIGENFUNCTIONS OF THE LATTICE 
HYPERBOLIC RUIJSENAARS-SCHNEIDER SYSTEM WITH 
EXPONENTIAL MORSE TERM 

J.F. VAN DIEJEN AND E. EMSIZ 


Abstract. We place the hyperbolic quantum Ruijsenaars-Schneider system 
with an exponential Morse term on a lattice and diagonalize the resulting n- 
particle model by means of multivariate continuous dual q-Hahn polynomials 
that arise as a parameter reduction of the Macdonald-Koornwinder polynomi¬ 
als. This allows to compute the n-particle scattering operator, to identify the 
bispectral dual system, and to confirm the quantum integrability in a Hilbert 
space set-up. 


1. Introduction 

It is well-known that the hyperbolic Calogero-Moser n-particle system on the line 
can be placed in an exponential Morse potential without spoiling the integrability 
[n ILW) . An extension of Manin’s Painleve-Calogero correspondence links the par¬ 
ticle model in question to a multicomponent Painleve III equation m- Just as for 
the conventional Calogero-Moser system without Morse potential, the integrability 
is preserved upon quantization and the corresponding spectral problem gives rise 
to a rich theory of remarkable novel hypergeometric functions in several variables 
[IMIIQIEIIHL]. 

An integrable Ruijsenaars-Schneider type (g-jdeformation [RSI IRl] of the hy¬ 
perbolic Calogero-Moser system with Morse potential was introduced in [S] and in 
a more general form in |D21 Sec. II.B]. Recently, it was pointed out that particle 
systems of this kind can be recovered from the Heisenberg double of SU{n,n) via 
Hamiltonian reduction [^. In the present work we address the eigenvalue prob¬ 
lem for a quantization of the latter hyperbolic Ruijsenaars-Schneider system with 
Morse term. Specifically, it is shown that the eigenfunctions are given by multi¬ 
variate continuous dual q-Hahn polynomials that arise as a parameter reduction of 
the Macdonald-Koornwinder polynomials [Kl IM2j . As immediate by-products, one 
reads-off the n-particle scattering operator and the commuting quantum integrals 
of a bispectral dual system [DGllG]. 

The material is organized is as follows. In Section [2] we place the hyperbolic 
Ruijsenaars-Schneider system with Morse term from [D2] on a lattice. The diago- 
nalization of the resulting quantum model in terms of multivariate continuous dual 
g-Hahn polynomials is carried out in Section |3l In Sections |4] and [5] the n-particle 
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scattering operator and the bispectral dual integrable system are exhibited. Finally, 
the quantum integrability of both the hyperbolic Ruijsenaars-Schneider system with 
Morse term on the lattice and its bispectral dual system are addressed in Section 


2. Hyperbolic Ruijsenaars-Schneider system with Morse term 

The hyperbolic quantum Ruijsenaars-Schneider system on the lattice was briefly 
introduced in [R2[ Sec. 3C2] and studied in detail from the point of view of its 
scattering behavior in [R4] (see also m Sec. 6] for a further generalization in terms 
of root systems). In this section we formulate a corresponding lattice version of the 
hyperbolic quantum Ruijsenaars-Schneider system with Morse term introduced in 
[D2l Sec. II.B]. 


2.1. Hamiltonian. The Hamiltonian of our n-particle model is given by the formal 
difference operator [D21 Eqs. (2.25), (2.26)]: 


E ( ( n 1 _ q.,-.. ) (2.1) 


H 


i=i 


-I- W-{ 


l<k<r. 


where 
w+{x) := 


(1 - <ig"^)(l - t 2 q^), w-{x) := J -^^(1 - 


tlt2 


qtots 


and Tj (j = 1,..., n) acts on functions / : K." ^ C by a unit translation of the jth 
position variable 


{Tjf){xi,. ..,Xn) = f(xi, . . .,Xj-i,Xj -f l,Xj+i,. . .,Xn). 

Here q denotes a real-valued scale parameter, t plays the role of the coupling pa¬ 
rameter for the Ruijsenaars-Schneider inter-particle interaction, and the param¬ 
eters tr (r = 0,..., 3) are coupling parameters governing the exponential Morse 
interaction. Upon setting tg = et^~^q~^ and U = e for r = 1,2,3, one has 
that w±{xj) —>■ when e —>■ 0. We thus recover in this limit the Hamil¬ 

tonian of the hyperbolic quantum Ruijsenaars-Schneider system given in terms 
of Ruijsenaars-Macdonald difference operators [Rl[ IMl) . By a translation of the 
center-of-mass of the form q^^ —^ cq^^ {j = 1,... ,n) for some suitable constant c, 
it is possible to normalize one of the t^-parameters to unit value; from now on it 
will therefore always be assumed that ts = 1 unless explicitly stated otherwise. 


2.2. Restriction to lattice functions. Let p-|-A:={p-|-A|AG A}, where A 
denotes the cone of integer partitions A = (Ai,...,A„) with weakly decreasingly 
ordered parts Ai > • • • > A„ > 0, and p = (pi,... ,p„) with 

Pj = {n- j)logg{t) (j = l,...,n). 


( 2 . 2 ) 
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The action of H ()2.ip (with ^3 = 1) preserves the space of lattice functions / : 
p T A —y Cl 

iHf){p + X)= Y. y+(x){fip + X + e,)-fip + X)) (2.3) 

l<j<n 
A+e j e A 

l<j<n 
A — e j e A 


where ei,..., Cn denotes the standard basis of M" and 


n+(A) 


(A) 


tlt2 


n 

l<fe<n 


1 — q^j~^k ’ 






) n 

l<k<r 


1 — t^~3 q^j~^k 


Indeed, given A S A, one has that v'^ (A) = 0 if A + Cj ^ A due to a zero stemming 
from the factor when Aj_i = Xj and one has that vJ{X) = 0 

if A — Cj ^ A due to a zero stemming from either the factor t — tq^^-^j+i when 
Xj = Aj+i or from the factor (1 — when A„ = 0. 


3. Spectrum and eigenfunctions 

Ruijsenaars’ starting point in |R4j is the fact that the hyperbolic quantum 
Ruijsenaars-Schneider system on the lattice is diagonalized by the celebrated Mac¬ 
donald polynomials |M1[ Ch.VI]. In this section we show that in the presence of 
the Morse interaction the role of the Macdonald eigenpolynomials is taken over by 
multivariate continuous dual g-Hahn eigenpolynomials that arise as a parameter 
reduction of the Macdonald-Koornwinder polynomials [K1 IM2] . 


3.1. Multivariate continuous dual g-Hahn polynomials. Continuous dual q- 
Hahn polynomials are a special limiting case of the Askey-Wilson polynomials in 
which one of the four Askey-Wilson parameters is set to vanish [KLSl Ch. 14.3]. 
The corresponding reduction of the Macdonald-Koornwinder multivariate Askey- 
Wilson polynomials [K] IM2j is governed by a weight function of the form 


A(0 := 


1 TT I 

(2-)%i,4Jno<.<2(iVc«o 


2 

l<j'<fc<n 





2 


(3.1) 


supported on the alcove 


A := {(^1, ^2, • ■ • , Cn) s K” I TT > > ^2 > ■ • ■ > Cn > 0}, (3.2) 

where {x)m ■= and {xi,...,xi)m ■= {xi)m ■ ■ ■ {xi)m refer to the 

g-Pochhammer symbols, and it is assumed that 


g, tG(0,l) and G (— 1,1) \ {0} (r = 0,1,2). 


(3.3) 
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More specifically, the multivariate continuous dual q-Hahn polynomials A € 

A are defined as the trigonometric polynomials of the form 


^>^(0 = X! (ca.m G C) 


such that 


and 


ca.a = 


n 


^<X 




n 




Afc 




Ja 

Here we have employed the dominance partial order 

V/r, A S A : ^ < A iff E < E Xj for k = 1,.... 

and the symmetric monomials 

mxiO ■■= E A e A, 

vewx 


(3.4a) 

(3.4b) 

(3.4c) 

(3.5) 

(3.6) 


associated with the hyperoctahedral group W = Sn^ {1, —1}" of signed permuta¬ 
tions. 

The present choice of the leading coefficient Ca,a in Eq- (13.4b() normalizes the 
polynomials in question such that P\(ip) = 1, where p = (pi,... ,/5„) is given by 
Pj = (n —j) log(t)-|-log(to), j = 1, • ■ •, (cf. [D31 Sec. 6]. |M21 Ch. 5.3]). With this 
normalization, the orthogonality relations obtained as the degeneration of those for 
the Macdonald-Koornwinder polynomials [K] Sec. 5], [D3I Sec. 7], |M21 Ch. 5.3] 
read: 


where 


/ Px{i)p^.{s,)Hm 

Ja 


ifA = p, 

0 otherwise, 


Aa :=Ao 

l<j<n 


{ioiit"' ^)xj 


(3.7a) 


(3.7b) 


X 


n 

l<j'<fe<n 


l_t^-3qXj-Xk (tl + fc-i)A._A, 
l-tk-3 


and 


Aq := 


n( 

l<_j<n 


iq,P)oo 

(t)oo 


(irist^ '’)oo)- 

0<r<s<2 


(3.7c) 


3.2. Diagonalization. Let £^(p-|-A, A) denote the Hilbert space of lattice func¬ 
tions f : p + A ^ <C determined by the inner product 

(/)P)a := ^ /(p + A) 5 (p -I- A) a a {f, g € £^{p + X, A)), 

AeA 


(3.8) 
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with p and Aa as in Eqs. (12.21) and (I3.7all - (|3.7c|) . and let L^(A, A(^)d^) be the 
Hilbert space of functions f : A ^ C determined by the inner product 

{f,g)A--= [ fiOmMm {f,g {A, AiOdO), (3.9) 

Jk 


with A taken from Eq. dan. We denote by : p + A —>• C the lattice wave 
function given by 


^^{p + X):=P^{0 (CgA, AeA). 

(3.10) 

Then the orthoffonalitv relations in Eos. (I3.7all~(l3.7cl) imolv that the associated 
Fourier transform F : £^(p -|- A, A) —>■ A^(A, Ad^) of the form 

(Ffm := if, = Yf^P + ^'^^dP + A)Aa 

AgA 

(3.11a) 

(/ G £^(p-|-A, A)) constitutes a Hilbert space isomorphism with an inversion formula 
given by 

(F-V)(P + A) = (/,V'(P + A))^= / f{OMP+^)^iOd^ 

Jk 

(/GL2 (A,AdO). 

(3.11b) 

Theorem 1. Let E denote the bounded real multiplication operator 
L2(A,Ade) by {Efm -Mono with 

acting on f G 

Mo ■■= Y (2 cos(e,) - r-do - t^-nn ). 

l<_ 7 <n 

For 

(3.12a) 

to = g ti = iod, t2 = tod 

(3.12b) 


with q, t and tr in the parameter domain (1331) . the hyperbolic lattice Ruijsenaars- 
Schneider Hamiltonian with Morse interaction H (lA3l) constitutes a bounded self- 
adjoint operator in the Hilbert space £^(^ 3 -A, A) diagonalized by the Fourier trans¬ 
form F (j3.11a|l . (|3.11bp ; 

H = F-^ oEoF. (3.12c) 

Proof. It suffices to verify that the Fourier kernel (13.1011 satishes the eigenvalue 
equation Hip^ = i5(^)'0j, or more explicitly that: 

t+(A)(V'5(p +A + Cj) - V'c(p +A)) + 

l<j<n 
A+fij G A 

Y ^7Wi'^dp + ^ - - '^dp + ^)) = Ei07dp + 

l<j<n 

X-ej£A 

This eigenvalue equation amounts to the continuous dual g-Hahn reduction of the 
Fieri recurrence formula for the Macdonald-Koornwinder polynomials correspond¬ 
ing to Eqs. (6.4), (6.5) and Section 6.1 of |D3j . □ 

It is immediate from Theorem[I]that the hyperbolic lattice Ruijsenaars-Schneider 
Hamiltonian with Morse interaction H (12.311 has purely absolutely continuous spec¬ 
trum in £^(p -I- A, A), with the wave functions '(/’{, ■C S A in Eq. (I3.10|l constituting 
an orthogonal basis of (generalized) eigenfunctions. 
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Remark 2. For ^2 —>■ 0 the lattice Hamiltonian H (j3.12c|l becomes of the form 

^ 1 _ q^3 


j=l \ l<k<n ^ 

+ toil - ) ( n 1 ^ - £0, 


(3.13a) 


l<fc<n 

k¥=j 


with X = p + X and 


(3.13b) 


£0 :=^(toi"-' +to-'i'“”)- 

i=i 

Indeed, this readily follows from Eqs. (12.ip . (l3.12bP with the aid of the elementary 
polynomial identity (cf. Example 2. (a) of [Mil Ch. VI.3]) 


Ed-i) n 


l<fc<r; 

k^3 


J2izj+t--n- 

1=1 


4. Scattering 


1=1 


In this section we rely on results from |D4| . permitting to describe briefly how 
the n-particle scattering operator for the hyperbolic quantum Ruijsenaars-Schneider 
system on the lattice computed by Ruijsenaars [R4) gets modified due to the pres¬ 
ence of the external Morse interactions. Specifically, the scattering process of the 
present model with Morse terms turns out to be governed by an n-particle scattering 
matrix S{^) that factorizes in two-particle and one-particle matrices: 


*5(0:= H ^i^j - ^k)si^j + ^k) n '®o(^l)> 


with 


3ix) := 


(ge^Ete-^")c 

(ge-“,te“)c 


and soix) := 


(ge"“)o, 

(ge-2“)c 


n 

0<r<2 


{ir 


(tr-e“)c 


(4.1a) 


(4.1b) 


which compares to Ruijsenaars’ scattering matrix ni<2<fc<rt '^(0 ~ 
responding model without Morse interactions [R4] . 

To substantiate further some additional notation is needed. Let us denote by 
T-Lq the self-adjoint discrete Laplacian in i'^iA) of the form 

(Ho/)(A):= ^ /(A + e,)+ ^ /(A-e,) (/e ^^(A)), 


\-\~Gj GA. 


l<j<n 
X — e 1 ^ A 


and let 


n:= A^/'^iH + eo)A-^^'^, (4.2) 

with H and £o taken from (12.31) and (I3.13bl) . respectively. Here the operator : 
£'^ip + A, A) i'^iA) refers to the Hilbert space isomorphism 

(A1/2/)(A) := A\/^fip + A) (/ G fip + A, A)) 


(4.3) 
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(with A := ^). Then (by Theorem[T|) 

+ with := A^^Va^^/^ (4.4) 

>v 1/2 

where A : L^(A, Ad^) —5> L^(A) denotes the Hilbert space isomorphism 

(a'^VKO := Ai/2(0/(e) (/ G L^(A, AdO) (4.5) 

(and E ()3.12al) is now regarded as a self-adjoint bounded multiplication operator 
in L^(A)). Furthermore, one has that 

’^0 = ^0 ^{E + eo)^o, 


where : £^(A) ^ A^(A) denotes the Fourier isomorphism recovered from E in 
the limit q,t ^ 0, E ^ 0 {r = 0,1,2). Specifically, this amounts to the Fourier 
transform 


(^o/)(0 = E/Wa?W 

(4.6a) 

AeA 


(/ G £'^{A)) with the inversion formula 


= f mxdm 

Ja 

(4.6b) 


(/ G L^(A)) associated with the anti-invariant Fourier kernel 


^ ^ w^W 

where sign(ri;) = ei ■ • • ejiSign(tT) for w = (cr, e) G IT = S'„ k {1,-1}" and po = 
(n,n- 1 ,..., 2 , 1 ). 

Let C'o(Areg) be the dense subspace of L?{K) consisting of smooth test functions 
with compact support in the open dense subset Aj-eg C A on which the components 
of the gradient 

VA(C) = (-2 sin(ei),..., -2 sin(^„)), ^ G A 

do not vanish and are all distinct in absolute value. We define the following unitary 
multiplication operator S : L^(A, d^) ^ L^(A, d^) via its restriction to C'o(Ai.eg): 

[Sfm := S{w^0f{0 (/ G Co(Areg), (4.7) 

where wj G IT for ^ G A^eg is the signed permutation such that the components of 
w^\/E{^) are all positive and reordered from large to small. 

Theorem 4.2 and Corollary 4.3 of Ref. [D4) now provide explicit formulas for the 
wave operators and scattering operator comparing the large-times asymptotics of 
the interacting particle dynamics relative to the Laplacian’s reference dynamics 
giWo* continuous dual q-Hahn reduction of [D4l Thm. 6.7]. 

Theorem 3 (Wave and Scattering Operators). The operator limits 

—s- lim 

converge in the strong £‘^(A)-norm topology and the corresponding wave operators 
intertwining the interacting dynamics with the discrete Laplacian’s dynam¬ 
ics are given by unitary operators in £'^{A) of the form 

oEo, 
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where the branches of the square roots are to he chosen such that 


f n 1/2 ^ (ge^^)oo l(te'^)oo| 
|(ge“)oo| (ie“)oo 


and 


soixf^'^ = 


(ge^^^)oo -r-r |(fre^^)oo| 

l(9e2“)oo| (i;e“)oo 


The scattering operator relating the large-times asymptotics of for t —>■ —oo 
and t —>■ +00 is thus giuen by the unitary operator 

s := oSoTq. 


5. Bispectral dual system 

The bispectral dual in the sense of Duistermaat and Griinbaum [DGIIG] of the 
hyperbolic quantum Ruijsenaars-Schneider system on the lattice is given by the 
trigonometric Ruijsenaars-Macdonald g-difference operators [Rll IHU- This bis¬ 
pectral duality is a quantum manifestation of the duality between the classical 
Ruijsenaars-Schneider systems with hyperbolic/trigonometric dependence on the 
position/momentum variables and vice versa [R3) . which (at the classical level) 
states that the respective action-angle transforms linearizing the two systems un¬ 
der consideration are inverses of each other. As a degeneration of the Macdonald- 
Koornwinder q-difFerence operator m Eq. (5.4)], we immediately arrive at a bis¬ 
pectral dual Hamiltonian for our hyperbolic quantum Ruijsenaars-Schneider system 
with Morse term. 

Indeed, the continuous dual q-Hahn reduction of the q-difference equation satis¬ 
fied by the Macdonald-Koornwinder polynomials m Thm. 5.4] reads 

n 

HPx = ExPx with Ea -1) (AeA), (5.1a) 

i=i 

where 

n 

H = - 1) + Vj{-OT-,^ - 1)) , (5.1b) 

and 

no<r<2(l 

“ (l_e2*G)(l-qe2*?j) H i _ 

^ ^ ' l<fc<n 

Here Tj^q acts on trigonometric (Laurent) polynomials p(e®^L ..., e®^") by a q-shift 
of the qth variable: 

(Tj.gP)(e®‘^®,..., e®«") := p(e®fi,..., , qe®«^, ,..., e®«"). 

In other words, the bispectral dual Hamiltonian H (I5.1bl) . (l5.1cl) constitutes 
a nonnegative unbounded self-adjoint operator with purely discrete spectrum in 
L^(A, Ad^) that is diagonalized by the (inverse) Fourier transform F (l3.11aL 
(l3llbll : 

H = FoEoF-^, (5.2) 

where E denotes the self-adjoint multiplication operator in ^^(p-|-A, A) of the form 
{Ef){p A) := Exf{p-\-X) (for AeA and / e £^(p -I- A, A) with {Ef, Ef)^ < oo). 
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6. Quantum integrability 

In this final section we provide explicit formulas for a complete system of com¬ 
muting quantum integrals for the hyperbolic quantum Ruijsenaars-Schneider Hamil¬ 
tonian with Morse term on the lattice H (12.31) and for its bispectral dual Hamilton¬ 
ian H (I5.1b|) . (I5.1cll . This confirms the quantum integrability of both Hamiltonians 
in the present Hilbert space set-up. 


6.1. Hamiltonian. The quantum integrals for the hyperbolic Ruijsenaars-Schneider 
Hamiltonian with Morse term are given by commuting difference operators Hi,, Hn 
that are defined via their action on f € i^{p + A, A) (cf. |D21 Eqs. (2.20)-(2.23)]): 

{Hf){p + X):= (6.1) 

,z-|j+|-|j_|(^)K 7 +.j_ (A)/(p -I- a -I- ej_^_ - ej_) 

J+, C{l,...,n} 

j+nj_=0. |J+|+|J_|<z 
\+ej^ —ej_ £A 


(A G A, Z = 1,..., n), where ej := for J C {1,..., n}, := {1,..., n} \ J 

and 


Vj+ j^(A) = i-^IJ+l(l-^+l-i)+^|J-l(|J-|-i) 
j&J+ * 


n \/^(i - 


ieJ- 


n(- 

i6J+ 

keJ- 




-){ 


j&J+ 

k^J+UJ- 


1 — J \ X — fk—jq\j—Xk+i J 

t — 

X — q^j ~^k ’ 


J--L ^ — -^k 3 Q^j 


n 


j&j- 

k^J+UJ- 


Uk.pW = 

^ f 

I+J-CK 

/+n/_=0.|/+|+|/_|=p 


Vi 6 A 


n \/^(l - 9^0(1 - 


At V 9^0 


nc 


_ f^~^k—j^Xj—Xk \ / ]^ f^ — ^-\-k—j^Xj—Xk-\-i- 


kGl- 


l — tk-jq\j-\k J\ I — fk-jq\j-\k+i ) 


-)( 


n 




TT - 

X — t^~3 Q^j ^k J. J. X — t^~3 Q^j ^k 

3^1+ ^ ^ 

k€K\(I+Ul-) feeA\(7+U/_) 


Xi-\k 
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For I = 1 the action of Hi (16.ip is seen to reduce to that of H (12.31) . The diagonal- 
ization in Theorem [T] generalizes to these higher commuting quantum integrals as 
follows. 


Theorem 4. For parameters of the form as in Theorem{ll the difference operators 
Hu...,Hr, dSTD constitute bounded commuting self-adjoint operators in the Hilbert 
space £^(/9 +A, A) that are simultaneously diagonalized by the Fourier transform F 

(IrTTaP . (IrlTbll ; 

Hi=F-^oEioF {l = l,...,n), (6.2a) 

where Ei denotes the bounded real multiplication operator acting on f G T^(A, Ad^) 

by {Ejm ■■= Eiiom with 

EliO := (6.2b) 

^ (2 cos(e,J - • • • (2 cos(C,J - H-% - ')• 

Proof. The eigenvalue equation Hiip^ = reads explicitly 

UjinJ^,i-\J+\-\J-\{^Wj+,J-W'^i{p + A + ej^ - ej_) 

J+,J_C{l....,n} 
j+nj_=0. |J+|+|J_|<i 

\+ej^ —ej_ GA 

= Ei{£.)i’dP + ^)- 

This eigenvalue identity corresponds to the continuous dual q-Hahn reduction of the 
Fieri recurrence formula for the Macdonald-Koornwinder polynomials in [D31 Thm. 
6.1], where we have expressed the eigenvalues Ei(f) in a compact form stemming 
from [KNSi Eq. (5.1)] (cf. also |DE[ Sec. 2.2]). □ 


6.2. Bispectral dual Hamiltonian. The continuous dual q-Hahn reduction of 
the system of higher g-difference equations for the Macdonald-Koornwinder poly¬ 
nomials in |D3I Sec. 5.1] reads 

HiPx = Ex,iPx (A e A, / = 1,..., n), (6.3a) 

where 

Exj := - e-E)... 

(6.3b) 

(cf. [Ki^ Eq. (5.1)]), and 

Hi := Y Ujui-iJiVejtj,<i, (6.3c) 

, 7 C{l.....rt}. 0 <|J|<i 
ei6{l -1}J67 


with Tej,q:= and 




n 


1 _ I 


jGJ '' yv 1 / Jgj 

k^J 


n 

3,k^J 

j<k 


X _ edej?i-l-e(=?fc) X _ qgi{f-j^j+i^kS,k) ’ 


— ^k) 
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Uk,p = (- 1 )^ 


n 


n 


0<r<2 


(1 - ) 


I<ZK,\I\=p ''jel 
€jC{l, —l},jG/ 


n 


(1 — 


jei 

keK\I 




Ht 

j,k^I 

j<k 


1 _ ie*(eiSj+«»!Sfc) I — +^3ik) 

— e®(<=j53+'=fc^fc) 1 — qe^i<^jii+'^kik) 


For 1 = 1, this reproduces the continuous dual q-Hahn reduction of the Macdonald- 
Koornwinder q-difference equation in Eqs. (I5.1al) - (l5.1cp . 

The q-difference operators Hi ,..., -ff„ extend the bispectral dual Hamiltonian 
H (I5.1bl) - (l5.1c[) into a complete system of commuting quantum integrals that are 
simultaneously diagonalized by the multivariate continuous dual q-Hahn polynomi¬ 
als. 


Theorem 5. For parameter values in the domain (USD, the q-difference operators 
Fdi,..., Hn constitute nonnegative unbounded self-adjoint operators with purely dis¬ 
crete spectra in L^(A, Ad^) that are simultaneously diagonalized by the (inverse) 
Fourier transform F (Id.llal) . (I3.11b|) .' 

Hi = FoEioF-\ l = l,...,n, (6.4) 

where Ei denotes the self-adjoint multiplication operator in F A, A) given 
by {Eif){p F A) := Exjf{p F A) (on the domain of f € H{p -I- A, A) such that 
{Elf, Elf)A < oo). 

Notice in this connection that although the domain of the unbounded operator 
Hi in A^(A, Ad^) depends on I, the resolvent operators {Hi — zi)~^,..., {Hn — Zn)~^ 
(with 01 ,..., € C \ [0, +oo)) commute as bounded operators on H{K, Ad^), and 

the q-difference operators Hi ,..., Fin moreover commute themselves on the joint 
polynomial eigenbasis P\, A G A. 

Remark 6. To infer that the eigenvalues E\^i (I6.3bl) are nonnegative—thus indeed 
giving rise to a nonnegative operator Hi in Theorem [S] —it is helpful to note that 
these can be rewritten as (cf. |D3[ Sec. 5.1]) 

Ex,l = , ^q-^^ . . . , ; t^-\t\ t”"!) 

with 

Epn{z^l, ■ ■ ■ , ^n', yi, • ■ • , yn) ■— ^ ^ ( 1) ^k{z^l, ■ ■ ■ , Zln)hl — k{yi, ■ ■ ■ , yn)- 

0<k<l 

Here ek{zi,..., Zn) and hk{yi, ■ ■ ■, yn) refer to the elementary and the complete 
symmetric functions of degree k (cf. mi Ch. 1.2]), with the convention that 
eo = ho = 1. The nonnegativity of the eigenvalues now readily follows inductively 
in the particle number n by means of the recurrence (cf. m Lem. B.2]) 

Ei,n{q ~^", tq~^\ ..., ; t’'~\t^,= 

{q-^i - t^-^)Ei.i,n-i{tq-^^ r-i) 

+ Ei^n-i{tq ~^^, ■ • ■, ;t\... 
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and the homogeneity 

(^-^1 ; • ■ • ; tyi T • ■ • : t -^/,n (^1 ; • ■ • ; yii • • • •> yn) • 

Remark 7. The hyperbolic Ruijsenaars-Schneider Hamiltonian with Morse term 
(EH) can be retrieved as a limit of the Macdonald-Koornwinder ^-difference op¬ 
erator m- In this limit the center-of-mass is sent to infinity, which causes the 
hyperoctahedral symmetry of the Macdonald-Koornwinder operator to be broken: 
while the permutation-symmetry still persists the parity-symmetry is no longer 
present. Indeed, the limit in question restores the translational-invariance of the 
interparticle pair-interactions enjoyed by the original Ruijsenaars-Schneider model 
and gives moreover rise to additional Morse terms that are not parity-invariant. It 
turns out that most of our results above can in fact be lifted to the Macdonald- 
Koornwinder level, even though such a generalization is presumably somewhat less 
relevant from a physical point of view. Specifically, the scattering of the cor¬ 
responding quantum lattice model associated with the full six-parameter family 
of Macdonald-Koornwinder polynomials was briefly discussed in |D41 Sec. 6.4], 
its commuting quantum integrals can be read-off from the Fieri formulas for the 
Macdonald-Koornwinder polynomials in [D31 Thm. 6.1], and the pertinent bis- 
pectral dual Hamiltonian and its commuting quantum integrals are given by the 
Macdonald-Koornwinder (/-difference operator m and its higher-order commuting 
(/-difference operators [D31 Thm. 5.1]. 
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